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Abstract 

C^ . The impact of the dispersion of the transport coefficients on the structure of the energy distri- 

bution function for charge carriers far from equihbrium has been investigated in effective medium 

fi I approximation for model density of states. The investigations shows that two regimes can be 

■ observed in energy relaxation processes. Below a characteristic temperature the structure of the 

energy distribution function is determined by the dispersion of the transport coefficients. Thermal 

^D I energy diffusion is irrelevant in this regime. Above the characteristic temperature the structure of 

(N 

O 

o 



the energy distribution function is determined by energy diffusion. The characteristic temperature 
depends on the degree of disorder and increases with increasing disorder. Explicit expressions for 
the energy distribution function in both regimes are derived for a constant and an exponential 
density of states. 



I. INTRODUCTION 

A characteristic feature of many disordered materials is the strong dependence of their 
transport properties on the frequency of the apphed external electric field for low frequen- 
cies. In that disordered materials differ from ordered materials, which typically only show a 
pronounced frequency dependence at high frequencies. Also the character of the impact of 
an alternating electric field on the transport coefficients in disordered materials is different 
from that in ordered ones. Whereas in ordered materials an increase of the frequency always 
reflects in a decrease of the conductivity in disordered materials an increase of the frequency 
leads to an increase of the real part of the conductivity. The sensibility of the transport 
coefficients to a change of the frequency increases with increasing amount of disorder. There- 
fore, a particular strong frequency dependence is observed in the strongly localized regime, 
in which transport proceeds by hopping between localized states, such as impurity bands, 
Anderson insulators, glasses, and polymers (see, e.g, Refs. l|y|3|y|5|). 

The same frequency dependence, which governs the response to an alternating electric 
field, determines also the relaxation properties of the system. Accordingly, the dispersion of 
the dynamical conductivity in the hopping regime reflects itself also in relaxation experiments 
close to equilibrium. It manifests itself, e.g., in transient current experimentsi*^, which test 
the diffusion properties of particle packets in a random environment. In such experiments 
the dispersion of the transport coefficients leads to dispersive transport, this is to deviations 
from the conventional relationship (_R^(t)) oc t for t < tperc, where (i?^) is the mean squared 
deviation of a particle packet at time t, which was a delta-pulse at time t = 0, and tp^rc is the 
percolation time^. Accordingly, the diffusion process breaks up into two regimes, a regime of 
anomalous diffusion, which takes place for t < tperc, and a regime of normal diffusion, which 
takes place for t ^ tp^rc- Since the percolation time depends exponentially on the disorder 
parameter it depends strongly on the material in question, i.e. there is no typical time scale 
for it in the experiments. In practice time scales between 10~^s and 5s have been observed 
(see discussion of experiments on p. 216 of Ref. [l|). 

Far from equilibrium mainly the effective rate analysis of Ref. |7| has been applied to 
analytical investigations of relaxation phenomena. This method assumes that the relaxation 
process far from equilibrium is not affected by the dispersion of the transport coefficients, 
and that therefore, the transport coefficients to be used in modeling the relaxation are 



the same as those at zero frequency. We would hke to mention that despite this fact also 
dispersive transport has been discussed in applications of the effective rate analysis, but that 
the nature of the dispersive transport discussed there, e.g. in Ref. p], is very different from 
that in transient current experiments close to equilibrium. Whereas the dispersive relaxation 
close to equilibrium results from the frequency dependence of the transport coefficients, and 
thus reflects the peculiarity that the charge carriers in hopping systems try to optimize 
their jumps, the dispersion in the effective rate analysis results solely from the fact that the 
particles move in an inhomogeneous distribution of sites in the four dimensional position- 
energy space. In this case dispersion arises since the transport coefficients depend on the 
instantaneous position of the particles in energy space, which changes with time. Therefore, 
similar dispersive contributions are also obtained for the relaxation of charge carriers in 
ordered systems. Clearly, they are much larger in the hopping regime than in ordered 
systems, since the transport coefficients depend much stronger on the energy in the hopping 
regime than in ordered systems. 

Given that the analytical results of the effective rate analysis correlate well with the 
results of numerical experiments for systems with strongly energy dependent density of 
states the question arises whether the dispersion of the transport coefficients is important in 
relaxation processes far from equilibrium or not, and if so, which quantities are determined 
by it. To answer this question at zero temperature we have derived an equation for the 



energv distribution function in the hopping regime in Ref. [Si] and investigated it further in the 
Refs. J2| and 1I| • The most important property of this equation is this the shape of the energy 
distribution function is entirely determined by the dispersion of the transport coefficients. 
The dispersion determines the first moment of the energy distribution function for systems 
with weakly energy dependent density of states and the second and the higher moments both 
for systems with weakly energy dependent density of states and for systems with strongly 
energy dependent density of states. These moments can be investigated experimentally by 
line width measurements in time resolved photoluminescence spectroscopy experiments. In 
an exponential density of states the dispersive nature of the energy relaxation leads to a 
soliton like motion of a particle packet in energy space at zero temperature. 

It is natural to ask whether these results apply also to systems at finite temperatures 
or not. If so, these results would provide a further characteristic, which would set disor- 
dered hopping systems apart from other physical systems. The latter generally share the 



feature that the structure of the energy distribution function is determined by energy diffu- 
sion processesiiii^. In the diffusion approximation, e.g. this fact reflects in that the second 
derivatives with respect to the energy are as important as the first derivatives and therefore 
can not be ignored^^. Below we investigate this question analytically for charge carriers far 
from the Fermi energy in a constant density of states and charge carriers in an exponential 
density of states. For these model densities of states we show that there is a characteristic 
temperature, which sets two regimes apart, a regime where the structure of the energy dis- 
tribution function is determined by the dispersion of the transport coefficients and a regime 
where the structure of the distribution function is governed by the energy diffusion current. 
The characteristic temperature discriminating these regimes is finite, depends on the amount 
of disorder, and increases with increasing disorder. 

II. THE MODEL 

We consider strongly localized non-equilibrium charge carriers far from the Fermi-energy. 
Such charge carriers can, e.g., be produced by illumination. After their excitation the charge 
carriers relax by phonon-assisted hops. We assume that the number of charge carriers is 
small, so that we can neglect Fermi-correlation. In this case the transport can be described 
by the simple rate equation 

Y^[nnWnm - riraWmn]- (l) 



dt 

Here n^ is the probability to find a particle on the site m with site energy e^ and position 
vector Rm, 

Wnm = d{uj - |e„m|)z/exp(-2a|i?„m| + -{enm - |enm|)), (2) 

^nm = ^n — ^m, Rnm = Rn — R-m, C(~^ IS the localization length, u is the attempt-to-escape 
frequency, and f3 is the inverse temperature. Both the positions of the sites and the site 
energies are random quantities. We assume that the sites are distributed randomly in space 
and that the site energies are distributed according to a distribution function, which is 
proportional to the density of states A^(e) of the system. 

In writing down Eq.Q we have also assumed that there is a maximal amount of energy 
u, which can be transferred from the electron to the phonon system, or vice versa, in one 
hop. The introduction of this energy scale is based on the observation that the electron 



phonon coupling constant tends rapidly to zero for phonons with large momentum q. The 
electron phonon coupling constant is a function which depends on both the overlap between 
the electron wave function and the phonon wave function and on the the Fourier transformed 
electromagnetic potential, which provides the coupling between the electron and the phonon 
wave. The overlap integral tends to zero for phonons with wave vector q > 2a, and thus 
renders the modes with q > 2a ineffective. The Fourier transformed potential is usually 
considered as a constant, in line with the deformation potential approximation. From the 
physical point of view this means that the electrons can move in such a way that they can 
keep the system neutral. In the strongly localized regime, however, there are no electrons 
which can move fast enough to screen out the alternations of the electric field produced by 
the phonon wave for most q. Accordingly, the Fourier-transformed Coulomb potential is a 
function which drops rapidly to zero already for tiny q. Due to these facts hops with large 
energy transfer can happen, but are not characteristic. In our approximation we therefore 
ignore such jumps. Since at present there are no further investigations in the literature on 
the energy scale u we take the point of view that a; is a phenomenological parameter, which 
can be determined by an experiment. 

In our investigation we focus on energy relaxation processes. Quantities of interest are, 
e.g., the mean energy 

(6)(t) = (5:6„n„(t)), (3) 

n 

and the mean squared deviation 

At) = {Y.^lMt))-{em (4) 

n 

of a particle packet in energy space. Here the bracket symbolizes the configuration average. 
Both quantities can be measured in time resolved photoluminescence experiments. In order 
to calculate these quantities we first assume that at time t = there is a particle packet, 
characterized by {^^(t = 0)}, and then use the function Pm'm (see Ref. |l[), the Green 
function of Eq.dH), to calculate the electron distribution for t > according to the equation 



n. 



,(t) = 5]n„(t = 0)P„„(t). (5) 



III. THE ENERGY DISTRIBUTION FUNCTION 

In order to derive an equation for the energy distribution function we use the relationship 

F(e', e|t) = f dRF{R\e', e, t) (6) 

between the energy distribution function F{e', e\t) and the diffusion propagator F{R\e', e, t). 
The latter is calculated from the configuration average of the Green function Pnm- To 
calculate the average we first introduce continuous coordinates according to the definition 

PiV; ^) = E ^i^' - Vn)PnmS{V - Vj , (7) 

nrn 

where V = (-R, e) and Vm = {Rm,^m)- Thereafter, we apply the effective medium- 
approximation by Gouchanour et. al.— . Although this technique has been developed origi- 
nally for investigations of hopping systems with topological disorder it can be checked that 
the extension of this method to systems with both topological and energetic disorder only 
a,„ou,.s to a change of notation. Accordingly, we can directly use the results of Ref Q. as 
discussed also in Ref.y^. 

n 

The effective medium approximation of Ref . [1^ reduces the calculation of the diffusion 
propagator to the solution to the system of integral equations 

{P{R,e'\0,e))=N{e')F{R\e',e), (8) 



sF{R\e',e)=6{R' -R)6{e' -e)+ fdpi [ F{\R - Ri\\e',ei)W{Ri\eue)N{e) 

- F{\R\\^',^)W{Ri\e',e,)N{ei)], (9) 

WiR\e' e) = ^(^'-'^) (10) 

^ ' ' ^ 1 + f{e)W{R\e, e') + f{e')W{R\e', e) ' ^'^^ 



and 



/-1(e) = s+ I deidRW{R\e, ei)N{ei). (11) 



Here s is the Laplace frequency, which corresponds to a Laplace transformation with respect 
to time, W{R\e', e) = Wnm\R„,n=R,e„=e' ,e^=e, and W is the renormalized transition probability. 
/(e) is the effective- medium constant, which for the present system is strictly speaking not 
a constant, but a function of energy. According to these equations the functions F and W 



depend also on s, but the dependence of F on s and the dependence of W^ on s have been 
suppressed to simphfy the notation. 

Despite the effective-medium approximation the calculation of the diffusion propagator 
is still a difficult task, since it requires finding solutions of integral equations. Due to the 
lack of symmetry, however, it is impossible to solve these integral equations analytically if 
the density of states depends on energy, as it is the case for most relevant physical systems. 
To simplify the equations we restrict the consideration to the limit jSuj <^ 1 and use the 
concept of quasi-elasticity. The introduction of this concept relies on the notion that the 
energy scale uj is small, so that 

•^ ^ ^ < 1. (12) 



/(e) de 
Then the renormalized transition probability takes the form 

W{R\e',e) = e{e' -e)e{uj~e' + e)W{R\e') 

+ e{e-e')e{uj-e + e')W{R\e)e-^^'-''\ (13) 

where 

and W{R) = z/exp(— 2ai?). If we furthermore restrict the consideration to high tempera- 
tures, where Pu <^ 1, and use the inequality (fT^ we find that the self-consistency equation 
flllll simplifies considerably. For s = we find that the solution of this equation is given by 

and for small 5, that is for s satisfying \pc{^, 0) — Pc{e, s)\ <C Pc{e, 0), we obtain the equation 
(Pc(e, 0) - pc(e, s)) exp(pc(e, 0) - pc{e, s)) = — --. (16) 

Here 

p,{e,s) =ln{2 f{e,s)u) (17) 

is the dimensionless characteristic hopping length, 

^o{e) = -^^exp{-p,{e,0)), (18) 

and Sd is the solid angle in d dimensions (5*2 = 27r, 5*3 = in). The dimensionless character- 
istic hopping length pc{e, s) and the characteristic hopping length Rc{e, s) are connected by 
the relationship pc{e, s) = 2aRc{e, s). 



In the same approximation the equation for the energy distribution function takes the 
form 

.f (,. ., . ,(, _ ,) + .||^,(,).„(, .)(|^ + P^)]. (19) 

Here 

2 

v{e,s) = ^A^(e) f dRW{R\e) (20) 

is the spectral energy relaxation rate. Note that, due to the Eqs. fll()|l and (jllll the effective 
transition probability depends on the Laplace frequency s. Therefore, the energy relax- 
ation rate is also frequency dependent, so that we can expect that dispersion of the energy 
relaxation rate manifests itself also in the moments of the energy distribution function. 

In order to find a concrete expression for the spectral energy relaxation rate we use 
Eq. ljlHj) and restrict the investigation to low frequencies. Doing so, we obtain 

v{e,s) = -uexp{-pc{e,s)). (21) 

The Eqs. (fT^ . (fTB|) . (fT^ and (pT|) yield a closed set of equations for the calculation of the 
energy distribution function, which can be used for the investigation of energy relaxation 
processes in an arbitrary density of states, provided the inequalities ()12j) and (3uj ^ 1 are 
satisfied. Below we apply this set of equations to the investigation of energy relaxation 
processes in model densities of states. 

IV. ENERGY RELAXATION IN A CONSTANT DENSITY OF STATES 

In this section we focus on energy relaxation processes in an unbounded constant density 
of states with Nuj sites per volume. Investigations of energy relaxation processes in such den- 
sities of states should be relevant for charge carriers in certain t ype s of Anderson-insulators 
above the glass temperature, as, e.g., for those used in the Refs.jl4l] and llq, for which field 
effect measurements indicate that the density of states varies only weakly with energy. 

Since for a constant density of states the transport coefficients are independent of energy 
the calculation simplifies considerably. In this case the differential equation (|T^ can be 
solved exactly. Its solution takes the form 



F{eo,e;s) 



exp(-f(e 



eoj 



4 "•" ujv{s) 



eo-e|) 



2u;f (s) 



§1 + « 

4 ajD(s) 



(22) 



To find the dependence of tlie distribution function on time we liave to perform an inverse 
Laplace transformation, according to 



F{eo,e-t) = ^Jjse''F{eo,e-s). 



(23) 



This integral yields the probability to find a particle on a site with site energy e at time t, 
if it was on a site with site energy eo at time t = 0. In calculating this integral we restrict 
the consideration to large times. For large times we expect that the particle packet has 
already moved down in energy space. It has somewhere in energy space a maximum and 
some width. The position of the maximum is characterized by the mean energy (e)(t) and 
the width by the dispersion <J^(t). We assume that t is so large, that eo > (e)(t) — <j{t). In 
this case we can use the saddle point approximation. The saddle point satisfies the equation 

d 



t 



P 



..^--^)(2 



P' 



+ 



(24) 



\ 4 ujv{s) 

The solution of this equation is a function so(eo,e, t). If we assume that the integral is 
determined by the saddle point we deduce that 



F(eo,e,t) ocexp((-- 



/3' 



+ 



-50(^0; ^,t) 



-)(eo-6)). 



^ 4 ujv{so{eo,e,t))' 
Accordingly, the maximum of the particle packet satisfies the equation 



(25) 



P' 



So{eo,e,t) 



de^^2 V 4 Luvisoieo,e,t)) 



)(eo-e)] 



-f- 



^ + 



soieo,e,t) 



(26) 



^ 4 ujv{so{eo,e,t))' 

Here we have used the fact that sq satisfies the Eq. fl24|) . Thus, we know that the characteristic 
So in the vicinity of the maximum is close to zero. Accordingly, in order to calculate the 
integral we can restrict the consideration to small s, and use this knowledge to expand the 
root in the exponent of Eq. (j22j) . Doing so, we obtain 

2s' 



t^ 



\ 4 ujv{s) 



^-il + ^ 

2^ P'ujviO) /32cut;(0)u;o 

2s' . 



{p'uj)'v'{0)' 



(27) 



Here it is important to realize that the third and fourth term of this expansion is of different 
physical origin. Both terms, as can be seen further below, contribute to the mean squared 
deviation. However, while the third term results entirely from dispersion the fourth term 
results entirely from thermal energy diffusion. Thus the relationship between these two 
terms determines, whether thermal activation affects the shape of the energy distribution 
function, or not. If we use the expansion (j^Tj) the integrals become Gaussian. Doing the 
integrals we find that 

F(eo, 6; t) = exp(-^^ J ^ ^ ^ ) 28 

where 

a\t)=2ujv{id)t{l + ^^^^^^). (29) 

Accordingly, at large times the particle packet moves with constant velocity. An increase of 
the temperature reflects in a decrease of the energy relaxation rate. The shape of the packet 
is Gaussian, and its mean squared deviation is given by Eq. ()29|) . In this equation the first 
term in the bracket on the right hand side results entirely from thermal diffusion and the 
second entirely from dispersion. Remarkably, the thermal contribution to the packets width 
is independent of temperature. Since 

its contribution is only important if the temperature is large, that is, if T ^ Tc, where 



kTc = ujJpc/{Qd). Accordingly, the second derivative in Eq. (fT!I|l is negligible for T <; Tc. 
Consequently, the relaxation is determined by the first order differential equation 

sF{eo, e) = 6{eo - e) + ^(i^(eo, e)uPv{s)), (30) 

for T <^ Tc. If the relaxation at zero temperature is also described percolation like then this 
equation differs from the equation for the calculation of the energy distribution function at 

only in that f (0) is replaced by ujf3v{0), so that all of the 
results derived in Ref . |^ are also valid in this regime, after a trivial change of parameters. 
Since Tc oc y/p^ the range of applicability of these results increases with increasing disorder. 



zero temperature (see Ref. 
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V. ENERGY RELAXATION IN AN EXPONENTIAL DENSITY OF STATES 

In this section we investigate the relaxation in an exponential density of states of the 
type 

iV(e) = iVoexp(|). (31) 

Here A^o and A are parameters characterizing the density of states. The consideration of the 
relaxation in such a density of states is relevant for investigations on relaxation phenomena 
in band tails of several amorphous semiconductors, as e.g., for the conduction band tail in 
amorphous hydrogenated silicon (see, e.g., Ref. J7[). 

Since the strong dependence of the density of states on energy manifests itself also in the 
transport coefficients it turns out to be hard to find an analytical solution for the differential 
equation ()19p in this case. Therefore, we focus at first on the solution for a particle which 
moves in a region with 

kT^P^ « 1. (32) 

That is, we first restrict the consideration to the shallow states in the tail. The lower limiting 

,.,.... ......„....,.._........,. »H,a- 

|l7l |. In the literature this region is also called the hopping down regime^. 

Within the range of the inequality ()32j) we can use the WKB- approximation (see the 
details of the calculation in appendix E)). Doing so, we find that the energy distribution 
function takes the form 

Fieo, e- 1) = . ^ exp(-^i^i|M)!) (33) 

at large times, where the mean energy is given by 

(e)(t)^-Aln(lln(|^^t)), (34) 

Po 3A 

(here po = Pci.^) exp(e/A)) and the mean squared deviation by 

^\t) = ^lermit)+^lspit)- (35) 

Here 

''''^™^^^ "^ Pc((e)(t))' ^^^^ 
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and 

Eq. p3|l differs from tlie predictions of tlie energy loss flopping model, which have been 
derived with the assumption that the relaxation at zero temperature is percolation like, 
in two points. First, while in the energy loss hopping model the energy relaxation rate is 
independent of temperature^ the energy relaxation rate is determined by the parameter 
^ijj^vtj IS. at finite temperatures in the limit (iuj <^ 1. Second, while the mean squared 
deviation of the particle packet in the energy loss hopping model is time independent, the 
dispersion becomes time dependent at finite temperatures. The time dependence results 
from the thermal contribution to the mean squared deviation. This contribution to the 
mean squared deviation is dominant if 

kT » cuv'p,((e)(t))/(6rf), (38) 

that is for high temperatures. In this temperature range the width of the energy distri- 
bution function increases with temperature. However, since the particle is sinking down 
its hopping length increases. From Eq. ljH^ we deduce that this increase tends to suppress 
the thermal contribution to the packet's width, and thus renders the disperse contribution 
more important. In contrast to the thermal contribution the disperse contribution is time 
independent, so that the packet moves in energy space without further distortion, once the 
thermal spreading of the packet has become insignificant. In this limit the motion of the 
packet in energy space is soliton like. Here an increase of the temperature results in a 
decrease of the packet width. 

The fact that the thermal spreading of the particle packet becomes unimportant outside 
the validity range of the inequality (|38p indicates that in this region the second derivative 
with respect to the energy in Eq. (fTT?|l is negligible, so that instead of Eq. (|T!Hl the simple 
equation 

sF(eo, e) = b{t - eo) + ^(^(^0, e)/5u;t;(e, s)). (39) 

can be used. 

In the derivation of our results we have used the inequality p2|) . which controls the validity 
of the WKB-approximation. However, the results of this calculation show that thermal 
spreading becomes unimportant when time goes by, that is when the hopping length becomes 
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large. From our point of view this is intuitively clear, since the impact of the dispersion 
of the transport coefficients in the hopping regime increases strongly with disorder. This 
fact suggests, that Eq. (jHU|) also holds outside the validity range of the inequality (|H^. so 
that Eq. (j39p can probably also be used for deep states in the tail, if only the condition (J12j) . 
which governs the range of applicability of the quasi-elastic approximation, is met. 

VI. CONCLUSIONS 

Our investigations show that the structure of the energy distribution function in the 
hopping regime is determined by the dispersion of the energy relaxation rate in a large 
temperature range. The characteristic temperature, governing the impact of the dispersion 
on the structure of the energy distribution function, depends on the strength of the disorder, 
and therefore on the details of the density of states. In both of the cases studied the impact 
of energy diffusion processes on the structure of the energy distribution function were only 
significant for kT ^ LO^J^Tj^e^ 0). This sets the situation in the hopping regime apart from 
that on the extended side of the metal-insulator transition, in which the structure of the 
energy distribution function is entirely determined by energy diffusio n^^i^^ . 

In a constant density of states the parameter pj^t^ 0) is independent of e. Therefore, 
a global characteristic temperature exists, which governs the infiuence of the dispersion. 
Although in this case the mean squared deviation increases linearly with time for large times 
the difference between the energy diffusion driven structuring of the energy diffusion function 
and a structuring due to the dispersion of the transport coefficients manifest themselves in 
the temperature dependence of the second moment of the energy distribution function. The 
latter decreases quadratically with increasing temperature at large temperatures. 

In an exponential density of states the situation is more intricate. In this case there is no 
globally defined characteristic temperature, since the criterion kT ^ ijy/j)^{e, 0) depends on 
the instantaneous position of the packet in the tail. Since the particles are sinking down their 
dimensionless hopping length pd^, 0) increases when time elapses. Accordingly, the thermal 
spreading of the energy distribution function, which is important in the initial phase of the 
relaxation in the shallow states of the tail, is getting unimportant when time goes by. Deep 
in the tail the structure of the energy distribution function is entirely determined by the 
dispersion of the energy relaxation rate. Energy diffusion processes are negligible in this 



13 



case. In this limit the dispersive nature of relaxation process leads to a soliton like motion 
of the particle packet in energy, i.e., the particle packet moves in energy space without 
distortion. 

We would like to mention that in both cases studied an increase of the temperature results 
in decrease of the energy relaxation rate. This result is in line with the results the numerical 



inergj 
■Mi 



simulations in Ref. [idl for a constant density of states, which have been performed in the 
same limit. It sets the situation in the limit of small energy transfer apart from that in 
conventional variable-range hopping systems, in which the relaxation rate usually speeds up 
with increasing temperature^. 

APPENDIX A: CALCULATION OF THE INVERSE LAPLACE TRANSFORM 
FOR THE EXPONENTIAL DENSITY OF STATES 

Within the range of the inequality (j32p we can use the WKB-approximation. However, 
in applying the WKB-approximation the following difficulty appears. To lowest order, that 
is within exponential accuracy, the WKB-approximation immediately predicts 



F(eo,e) ocexp(--(e-eo) - I / dei 



3^ s 

If we use the inequality fl32j) we can generate a systematic expansion of the preexponential 
factor of the diffusion propagator with respect to the parameter kTp'^{e), where the prime 
indicates the derivative with respect to energy. However, since the error estimate in the WKB 
approximation is of multiplicative type, non of the approximations can satisfy probability 
conservation exactly. Therefore, we only use the exponent, and determine the preexponential 
factor by normalization. 

In investigating the energy distribution function further we again restrict the consider- 
ation to large times, where eg ^ {^){t) ~ cr(t). In this case the calculation of the time 
dependence of the energy distribution function amounts to the calculation of the integral 



m = ^ fdsexpist- rde,L^ + —^-^)). (A2) 

Ziri Jc Je \ 4 ujv[ei,s) 2 

In order to calculate this integral we use the same method as for the constant density of 
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states. Doing so, we find that 



(t-Q 



2 



m = I ,, , ^^p(-(^^7^))- (A3) 



^^13 — T^' (A4) 



Here 

/■eo 

75o;t'(ei,0) 



and 

where the thermal contribution is given by 



D\t)=Dl^^.^{t) + Dl(e), (A5) 



DtHerM-^l ^3^2^2(,^^0)' ^"^^^ 



and the disperse contribution by 

Now we use the fact that our calculation of the inverse Laplace transformation is only correct 
in the vicinity of the peak. Since the peak is situated at t = t^ the maximum of the packet 
in energy space (e)(t) is situated at e = (e)(t), where 

t = tm)- (A8) 

If we use Eq. ()A8|) to calculate the mean energy at large times we obtain Eq. ()34|) . Since our 
calculation of the inverse Laplace transformation is only correct in the vicinity of (e) (t) we 
now replace e by (e)(t) in Eg. ()A3j) and then adjust the preexponential factor, as discussed 
above. Doing so, we obtain Eq. (jHH|) . where 



aL.™(t) = 2DtHerm{{em)P'uj\\{e){t),0) 
kTA 



and 



pciiemr 



(A9) 



4.P = 2D,,,,i{em)pWv\{e)it),0) ^ -^. (AlO) 
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